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Abstract. The well-known gyrokinetic problem regards the perturbative expansion related to the dynamics of a charged 
particle subject to fast gyration motion due to the presence of a strong magnetic field. Although a variety of approaches have 
been formulated in the past to this well known problem, surprisingly a purely canonical approach based on Lie transform 
methods is still missing. This paper aims to fill in this gap and provide at the same time new insight in Lie-transform 
approaches. 



INTRODUCTION: TRANSFORMATION APPROACH TO GYROKINETIC THEORY 

A great interest for the description of plasmas is still vivid in the scientific community. Plasmas enter problems related 
to several fields from astrophysics to fusion theory. A crucial and for some aspects still open theoretical problem is the 
gyrokinetic theory, which concerns the description of the dynamics for a charged point particle immersed in a suitably 
intense magnetic field. In particular, the "gyrokinetic problem" deals with the construction of appropriate perturbation 
theories for the particle equations of motion, subject to a variety of possible physical conditions. Historically, after 
initial pioneering work [1, 2, 3], and a variety of different perturbative schemes, a general formulation of gyrokinetic 
theory valid from a modern perspective is probably due to Littlejohn [4], based on Lie transform perturbation methods 
[5, 6, 7, 8]. For the sake of clarity these gyrokinetic approaches can be conveniently classified as follows (see also 
Fig.l): 

A) direct non-canonical transformation methods: in which non-canonical gyrokinetic variables are constructed by 
means of suitable one-step [1], or iterative, transformation schemes, such as a suitable averaging technique [9], a 
one-step gyrokinetic transformation [10], a non-canonical iterative scheme [11]. These methods are typically difficult 
(or even impossible) to be implemented at higher orders; 

B) canonical transformation method based on mixed-variable generating functions: this method, based on canonical 
perturbation theory, was first introduced by Gardner [2, 12] and later used by other authors [13]). This method 
requires, preliminarily, to represent the Hamiltonian in terms of suitable field-related canonical coordinates, i.e., 
coordinates depending on the the topology of the magnetic flux lines. This feature, added to the unsystematic 
character of canonical perturbation theory, makes its application to gyrokinetic theory difficult, a feature that becomes 
even more critical for higher-order perturbative calculations; 

C) non-canonical Lie-transform methods: these are based on the adoption of the non-canonical Lie-transform pertur- 
bative approach developed by Littlejohn [4]. The method is based on the use arbitrary non-canonical variables, which 
can be field-independent. This feature makes the application of the method very efficient and, due to the peculiar 
features for the perturbative scheme, it permits the systematic evaluation of higher-order perturbative terms. The 
method has been applied since to gyrokinetic theory by several authors [14, 15, 16, 17]; 

D) canonical Lie-transform methods applied to non-canonical variables: see for example [18]. Up to now this 
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FIGURE 1. The transformation approach to gyrokinetic theory: |~2~| : see Gardner[2],[12]; [J]: see Hahm et al. [18] and present 
theory; [4]: first obtained by Alfven [1]; [T|: see Littlejohn [4]. 



method has been adopted in gyrokinetic theory only using preliminar non-canonical variables, i.e., representing the 
Hamiltonian function in terms of suitable, non-canonical variables (similar to those adopted by Littlejohn). This 
method, although conceptually similar to the developed by Littlejohn, is more difficult to implement. 

All of these methods share some common features,in particular: 

- they may require the application of multiple transformations, in order to construct the gyrokinetic variables; 

- the application of perturbation methods requires typically the representation of the particle state in terms of suitable, 
generally non-canonical, state variables. This task may be, by itself, difficult since it may require the adoption of a 
preliminary perturbative expansion. 

An additional important issue is the construction of gyrokinetic canonical variables. The possibility of constructing 
canonical gyrokinetic variables has relied, up to now, on essentially two methods, i.e., either by adopting a purely 
canonical approach, like the one developed by Gardner [2, 12], or using the so-called "Darboux reduction algorithm", 
based on Darboux theorem [4]. The latter is obtained by a suitable combination of dynamical gauge and coordinate 
transformations, permitting the representation of the fundamental gyrokinetic canonical 1-form in terms of the canon- 
ical variables. The application of both methods is nontrivial, especially for higher order pertubative calculations. The 
second method, in particular, results inconvenient since it may require an additional perturbative sub-expansion for the 
explicit evaluation of gyrokinetic canonical variables. 

For these reasons a direct approach to gyrokinetic theory, based on the use of purely canonical variables and transfor- 
mations may result a viable alternative. Purpose of this work is to formulate a "purely" canonical Lie-transform theory 
and to explicitly evaluate the canonical Lie-generating function providing the canonical gyrokinetic transformation. 



LIE-TRASFORM PERTURBATION THEORY 



We review some basic aspects of perturbation theory for classical dynamical systems. Let us consider the state x of a 
dynamical system and its <f -dimensional phase-space M endowed with a vector field X. With respect to some variables 
x = |x ! } we assume that X has representation [6] 

dx' 

where £ is an ordering parameter. We treat all power series formally; convergence is of secondary concern to us. By 
hypothesis, the leading term Xo of (1) represents a solvable system, so that the integral curves of X are approximated 
by the known integral curves of Xo. The strategy of perturbation theory is to seek a coordinate transformation to a new 
set of variables {*'}, such that with respect to them the new equations of motion are simplified. Since (1) is solvable 
at the lowest order, the coordinate transformation is the identity at lowest order, namely 

x i =x i + 0{e) (2) 

The transformation is canonical if it preserves the fundamental Poisson brackets. It can be determined by means of 
generating functions, Lie generating function or mixed-variables generating functions, depending on the case. In the 
Lie transform method, one uses transformations T which are represented as exponentials of some vector field, or rather 
compositions of such transformations. To begin, let us consider a vector field G, which is associated with the system 
of ordinary differential equations 

dx' : , , 

^ =<?(*), (3) 

so that if x and x are initial and final points along an integral curve (3), separated by an elapsed parameter e, then 
x, = Tx. In the usual exponential representation for advance maps, we have 

r = exp(eG). (4) 

We will call G the generator of the transformation T. In Hamiltonian perturbation theory the transformation T is 
usually required to be a canonical transformation. Canonical transformations have the virtue that they preserve the 
form of Hamilton's equations of motion. Canonical transformation can be represented by mixed-variable generating 
function, as in the Poincare-Von Zeipel method or by means of Lie transform. In the latter method vector fields G are 
specified through the Hamilton's equations. Following a more conventional approach, we can write the (3) in terms of 
the transformed point 

i-M 0> 

The components of the above relation are just Hamilton's equations in Poisson bracket notation applied to the 
"Hamiltonian" (Lie generating function) co, with the parameter e the "time." Equation (5) therefore generates a 
canonical transformation for any e to a final state x whose components satisfy the Poisson bracket condition 

[Qi,Qj\ = \pi,Pj\=0 (6) 

[9i,Pj] = Stj- (7) 

To find the transformation T explicitly, we introduce the Lie operator L= [©,...]. Recalling that coordinate compo- 
nents of vector are subject to pull back transformation law, then one gets 

dT 

- = -TL (8) 



with the formal solution 



T = exp 



J £ L(e')de' 



(9) 



For any canonical transformation the new Hamiltonian H is related to the old one by 



H = T~ l H + T~ l de'T (e') " ' . (10) 



To obtain the perturbation series one can expand (0,L, T,H and H as power series in e 



CO 



M = £ e n M„ 



(11) 



n=0 



where M represents (0,L,T,H. From (8), equating like powers of e, we obtain a recursion relation for the 
(On,L n ,T n ,H n (n > 0) which with 7q = 1, gives T„ in terms of L n and o„ in all orders. 



THE CANONICAL LIE TRANSFORM APPROACH TO GYROKINETIC THEORY 



The customary approach based on Lie-transform methods and due to Littlejohn [4] adopts "hybrid" (i.e., non-canonical 
and non-Lagrangian) variables to represent particle state, i.e., of the form z = (y,<j>). There are several reasons, usually 
invoked for this choice. In the first place, the adoption of hybrid variables may be viewed, by some authors, as 
convenient for mathematical simplicity. However, the subsequent calculation of canonical variables (realized by means 
of Darboux theorem) may be awkward and give rise to ambiguities issues [13]. Other reasons may be related to the 
the ordering scheme to adopted in a canonical formulation: in fact, in gyrokinetic theory, the vector potential A in 
the canonical momentum must be regarded of order I /0(e) while keeping the linear momentum of zero order, i.e., 
p =mv+^A. As a consequence, in a perturbative theory p must be expanded retaining at the same time terms of order 
I /0(e) and 0(e°), a feature which may give rise to potential ambiguities. According to Littlejohn [4] this can be 
avoided by the adoption of suitable hybrid variables, which should permit to decouple at any order the calculations 
of the perturbations determined by means of suitable Lie-generators. However, a careful observation reveals that the 
same ambiguity (ordering mixing) is present also in his method. In fact, one finds that the first application of the 
non-canonical Lie-operator method, yielding the lowest order approximation for the variational fundamental 1-form, 
provides non-trivial contributions carried by the first order Lie-generators. Probably for this reason, his approach is 
usually adopted only for higher-order calculations where ordering mixing does not appear. 

In this paper we intend to point out that canonical gyrokinetic variables can be constructed, without ambiguities, 
directly in terms a a suitable canonical Lie-transform approach, by appropriate selection of the initial and final 
canonical states (see path 13 in the enclosed figure), i.e., respectively x = (q,p) and X' = (Q[,P[, typ ,Py p ,(j)' f^). 
The latter are, by construction, gyrokinetic, i.e., the corresponding Hamiltonian equations of motions are independent 
of the gyrophase angle <p . We want to show that the transformation x — > X' can be realized, in principle with arbitrary 
accuracy in e, by means of a canonical Lie transformation of the form: 



being © = o(X',e) the corresponding Lie generator. In order to achieve this result, we shall start demanding the 
following relation between the fundamental differential 1 -forms, i.e., the initial and the gyrokinetic Lagrangians, which 
can be shown to be of the form: 




(12) 



dtL(x, ^-x,t) = dt3?(X', —X',t) + dSi + dS 2 + dS 3 + dS 4 . 



dt ' ' y 'dt 



(13) 



Here Si, S2,S^,S4 are suitable dynamical gauges functions, i.e., 




(14) 




(15) 



53 = ep' - my, 



(16) 




(17) 



where m and Ze are respectively the mass, the electric charge of the particle and ep' = — £ w ^ b the Larmor radius. 
Moreover, w' is a vector in the plane orthogonal to the magnetic flux line, while b' = B(r',t)/B(r',t), Q.' = ^fis 



the Larmor frequency and finally primes denote quantities evaluated at the guiding center position r'. In particular, 
vq = VQ(\j/ p ,Q' ,Q[,t) To the leading order in e one can prove 



r = r' + ep', 
v = v' + w'. 

The remaining notation is standard. Thus, up to 0(e) terms, there results 

v' = M 'b' + v^, 
w'=w' (e'j cos <j>' + e' 2 sin(j)') , 



(18) 
(19) 



w' = ^Wtf, 

mc 
Ze 



(20) 



pi _ mc ,,/ 



where \' E = cE' x b'/B' is the electric drift velocity and evaluated at the guiding center position and jii' is the magnetic 
moment, both evaluated at the guiding center position. Here we have adopted the representation of the magnetic field 
by means of the curvilinear coordinates ( y/ p ', Q' , Q\ ) where y/ p ' ,Q' are the Clebsch potentials according to which the 
magnetic field reads B' = Vy/^ x VQ', whereas we have introduced the covariant representation for the electric drift 
velocity \' E = v' v V\j/ p + v'qVQ'. The gyrokinetic Hamiltonian JT(x',f ), defined by means of 



r(x't\-p> dQ Wp' dw 'p +P ' d ^' V>(X' d x't\ 



(21) 



reads 



Ze 
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Here ST is the kinetic energy term, whereas canonical momenta read 
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(22) 

(23) 

(24) 
(25) 



Let us consider, for instance, the equation for P' A . We notice that p' coincides with gj.!' the first order Lie generator of 
the transformation (18). Therefore, PL results: 



(26) 



p * = - 1 ' * + i • { I A ' + \t v ' a ' • & + \ I V ' A ' + my ) - »**£ ) -£ ) = - 1 2»**£ ) -£ 

where, neglecting contributions of higher orders, the first term on the r.h.s. has been evaluated at the effective position 
r. Thus, denoting Pa = |t • ^i(x, j|x,f ) = — ^ • p T , the equation can be cast in the following form 

P'a=Pa + e[p'a,co], (27) 

where CO is the phase function: 

(0 = mj dQ'v Q ( V ' p ,Q>, Q\ , t) + ^A' • g<, 1} + Igg> • ^ V'A' • gg> + g« • mv+m£l' J ■ g\ (28) 




(29) 



and similarly 




dr 



with 



d\j/ p 



(30) 



Therefore, it follows that CO is really the Lie generating function of the canonical gyrokinetic transformation x — > X'. 
The calculation of CO is the sought result. In terms of CO the purely canonical gyrokinetic approach is realized. The 
procedure can be extended to higher orders to develop a systematic perturbation theory. 



Work developed in the framework of the PRIN Research Program "Programma Cofin 2002: Metodi matematici delle 
teorie cinetiche"( MIUR Italian Ministry) and partially supported (for P.N.) by the National Group of Mathematical 
Physics of INdAM (Istituto Nazionale di Alta Matematica), (P.N) by the INFN (Istituto Nazionale di Fisica Nucleare), 
Trieste (Italy) and (M.T.) by the Consortium for Magnetofluid Dynamics, University of Trieste, Italy. 



1. H. Alfen, Cosmical Electrodynamics, Oxford University Press, Oxford 1950. 

2. C.S. Gardner, Phys. Rev. 115, 791 (1959). 

3. T.G. Northrop and E. Teller, Phys. Rev. 117, 215 (1960). 

4. R.G. Littlejohn, J. Math. Phys. 20, 2445 (1979). 

5. R.G. Littlejohn, Phys.Fluids 24, 1730 (1981). 

6. R.G. Littlejohn, J. Math. Phys. 23, 742 (1982). 

7. R.G. Littlejohn, J. Plasma Phys. 29, 111 (1983). 

8. J.R. Cary and R.G. Littlejohn, Ann. Phys. (N.Y.) 151, 1 (1983). 

9. A.I. Morozov and L.S. Solov'ev, in Reviews of Plasma Physics, Edited by Acad. M.A. Leontovich (Consultants Bureau, New 
York, 1966), Vol. 2, p. 201. 

10. LB. Bernstein and P.J. Catto, Phys.Fluids 28, 1342 (1985). 

11. B.WeyssowandR. Balescu, J. Plasma Phys. 35, 449 (1986). 

12. J. Berkowitz and C.S. Gardner, Commun. Pure Appl. Math., 12, 501 (1959). 

13. H. Weitzner, Phys. Plasmas, 2, 3595 (1995). 

14. D.H.E. Dubin, J.A. Krommes, C.Oberman and W.W.Lee, Phys.Fluids 11, 569 (1983). 

15. W. W. Lee, Phys. Fluids 26, 556 (1983). 

16. T.S. Hahm, Phys. Fluids 31, 2670 (1988). 

17. A.J. Brizard, Phys. of Plasmas 2, 459, (1995). 

18. T.S. Hahm, W.W. Lee and A. Brizard, Phys. Fluids 31, 1940 (1988). 



ACKNOWLEDGMENTS 



REFERENCES 



